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A convenient measure of a map or flow's chaotic action is the topological entropy. 
In many cases, the entropy has a homological origin: it is forced by the topology of 
the space. For example, in simple toral maps, the topological entropy is exactly equal 
to the growth induced by the map on the fundamental group of the torus. However, 
in many situations the numerically-computed topological entropy is greater than the 
bound implied by this action. We associate this gap between the bound and the 
true entropy with 'secondary folding': material lines undergo folding which is not 
homologically forced. We examine this phenomenon both for physical rod-stirring 
devices and toral linked twist maps, and show rigorously that for the latter secondary 
folds occur. 
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I. INTRODUCTION 

In many industrial applications, fluid is stirred by moving rods to achieve homogeneity. 
Such a simple system also serves as a testbed for ideas about mixing and transport bar- 
riers. In the past few years, a topological description of rod-stirring has emerged (see for 
example [IHI]). In this framework, the two-dimensional fluid lives in a surface with holes, 
with the holes corresponding to moving rods and fixed baffles. This topological description 
applies to all fluids, but is most useful for very viscous flows. 

A consequence of the topological description is that we can compute lower bounds on the 
topological entropy of the system. The topological entropy is, roughly speaking, the growth 
rate of material lines in a fluid The lower bound on the topological entropy is a 
consequence of continuity - material lines cannot cross the rods, so they must grow at least 
as fast as dictated by the motion of the rods. 

We can then ask about the sharpness of this lower bound: we measure the actual growth 
rate of material lines in the flow, and compare it to the lower bound. In many systems, 
the 'gap' between the lower bound and the actual value is quite small - on the order of a 
few percent of the value. In other cases, the gap is very large, so that the observed value 
of the entropy is much larger than the lower bound. The central question of this paper is: 
what accounts for the difference in these two cases? In both cases, much of the entropy 
is 'homological' and arises from obstructions in the domain - the rods and baffles. This 
homological entropy depends only weakly on the details of the physical system, and can be 
deduced without actually solving the dynamical fluid equations. However, in the cases with 
a large gap it seems clear that there is another, dynamical source of entropy, which can 
depend more strongly on system parameters, such as the size of the stirring rods and the 
fluid being considered. We propose that the so called dynamical entropy is predominately 
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due to secondary folding, where material lines exhibit folds that are not directly associated 
with a topological obstacle. 

When the difference between the bound and the measured entropy is small, we refer to 
those bounds as sharp. We will investigate the sharpness of the bound by considering simple 
three-rod stirring devices. The motion of the rods is limited to a sequence of clockwise or 
counter-clockwise exchanges with one of their neighbors, following circular paths [TJ. These 
types of motions map naturally to generators of the braid group [?]: a clockwise exchange 
of the ith and (i + l)th rod is denoted o~i, and a counter-clockwise exchange by a^ 1 . A 
sequence of af l is called a stirring protocol. Since the fluid we consider is highly viscous, the 
stirring protocols define a periodic motion of the fluid elements, obtained by solving Stokes' 
equation for incompressible flow. Figure [T] shows some examples of stirring protocols, as 
well as their action on typical material lines. 

We will illustrate the gap between the bound and measured value using two examples in 
Section]!!} Because the disk with three punctures (rods) has such an intimate connection with 
the torus, we will examine toral maps in Section |III[ We observe the same gap phenomenon 
as in the fluid case, but because the map is a much simpler system we are able to prove some 
explicit results about toral maps: in particular, we demonstrate the presence of secondary 
folds, or 'kinks.' Finally, in Section [TV] we summarize our results and discuss future work. 



II. SOME OBSERVATIONS 

As an illustration of the gap phenomenon in hydrodynamic flows, consider a stirring 
device with three rods that initially lie in a line. There are many ways to move the rods 
and stir the fluid, but we are primarily interested in the topological aspects of the motion. 
As discussed in the introduction, we focus on sequences of interchanges of neighboring rods, 
which we label by generators of the braid group [U [2]. Here, we consider two stirring 
protocols, given in terms of generators by aicr^ 1 and ciaf , where generators are read left to 
right. The rod motions are depicted in the insets in Figure [1} The exact details of how we 
move the rods make no difference to the topology of the system, though it will in general 
affect the measured growth rate of material lines. The speed of motion is irrelevant, since 
we are in the Stokes flow regime. In the simulations below, we always use circular paths and 
constant speed. 



A. Computing /ihom and /ifl ow 

Two types of topological entropy will be computed here. The first, called /ihom, is the 
entropy coming from the homology, i.e. the entropy due to the rod motion alone. This is the 
lower bound alluded to in the introduction: it is the growth rate of a hypothetical 'rubber 
band' caught on the rods, and is independent of both the type of fluid being stirred and the 
specifics of the rod motion. This entropy is computed directly from the braid describing the 
rod motion. The second entropy we compute, h,R OW , is obtained by directly measuring the 
growth rate of material lines in the flow. We have ^hom — ^flow, since material lines cannot 
(asymptotically) grow slower than the minimum rate dictated by the rods. We now discuss 
how to compute h hom and /ifl ow . 

Computing /ihom can be difficult when dealing with a large number of rods or long se- 
quences of generators, though rapid techniques have been developed [H |9j. However, for 
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FIG. 1. The effect of two different stirring protocols on the same initial material line. The insets 
show the motion of the rods. The images are after four and three periods, respectively. In (b) 
boxes highlight some instances of secondary folding. There is no secondary folding visible in (a). 



three rods the entropy is very easy to compute. In that case, the topological entropy is 
equal to logarithm of the spectral radius (largest magnitude over all eigenvalues) of the 
Burau representation matrix [TUHH] -0 For braids of the form c^ovT^ the Burau matrix is 

'i + ki e 
k i 



l u 2 



(1) 



If /ihom > 0, then, under repeated iterations, a rubber band caught on the rods will grow in 
length at an exponential rate. In this case, the braid is called pseudo-Anosov. If h^ om = 0, 
the rubber band will not grow exponentially, and the braid is said to be finite- order. This 
terminology comes from the Thurston-Nielsen theory of classification of surface homeomor- 
phisms [To| ITT] . There is a third case in the classification, called reducible, which is not 
relevant here. Mixers with pseudo-Anosov stirring protocols are usually good at mixing; 
the exponential stretching and folding of material lines leads to a growth in the interface 
between solutes, which in turn allows diffusion or chemical reaction to act more rapidly. 

Now for /ifl ow : in theory it is computed by taking the supremum of growth rates over all 
material lines in the fluid [6]. In practice, a typical material line will eventually grow at a 
maximal rate, as long as some part of it is in the ergodic component with the fastest growth 
rate. Thus, we can get a good measure of h,R OW for each protocol by computing the rate at 
which a typical material line stretches. We do this by solving numerically for the motion 
of the individual fluid elements making up a material line, taking care to interpolate new 
points as the distance between neighbors becomes too large. 



B. Stirring Protocol a±a 2 1 vs. <J\o\ 

We compare two stirring protocols, given in braid form by aia^ 1 and o~\o~\. Using the 
matrix (JIJ, we find that they both have /ihom = log((3 + \/5)/2) ~ 0.962. Despite this, the 



effect on material lines is quite different, as can be see by comparing Figures 1(a) and 1(b) 



1 The (reduced) Burau representation is a representation of the braid group on n strings in terms of matrices 
of dimension n — 1. The representation arises from an action on homology on the double-cover of the 
punctured disk. See [7J [TS] for more details. 
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FIG. 2. Growth of the length of a material line for the two stirring protocols. Since they are both 
pseudo-Anosov systems, the slope, i.e. the growth rate tends to a constant - namely the topological 
entropy. The dashed line has slope equal to the lower bound, fohom = 0.962. 



Notice that for the 0\(s 2 protocol (Fig. 1(a)) the material line forms very smooth and 



regular layers; the only large visible fol ds are due to wrapping around the rods. On the 
other hand, for the <J\(j\ protocol (Fig. 1(b)), the material line has extra folds that are 
not due to wrapping around a rod, as highlighted by the boxes. These are what we call 
'secondary folds.' 

From simulation^] we can measure /ifl ow for each stirring protocol. Figure [2] shows the 
measured length of a material line during several periods of the flow. The length is plotted 
on a log scale, and the (asymptotic) slope of the line is h^ ow . In this way, we measure 
lij ow = 0.991 for the Oxo^} protocol, and 7tfl ow = 1.61 for the <J\<j\ protocol. For cricr^ 1 , the 
rod entropy and the flow entropy agree well (roughly a 3% difference). However, for 01 erf, 
the flow entropy is much larger than the rod entropy (/ihom accounts for roughly 61% of 
/ijow). We hypothesize that the larger measured entropy is due to additional stretching of 
the material line due to the secondary folds. 



III. TORAL LINKED TWIST MAP 



It is problematic to define secondary folding rigorously in a practical context, which also 
makes it difficult to identify its causes. To make headway, we look at a simpler system 
- a map instead of a hydrodynamic flow - which can be analyzed more thoroughly. The 
class of maps we'll examine are toral linked twist maps. Linked twist maps (or LTMs) have 
been studied extensively [T81 - I25] . They are non-uniformly hyperbolic, and are especially 
relevant here because of an intimate connection to the three-rod stirring protocols above. 
The connection arises through the orientable double cover [H [26] : the torus can be regarded 
as a double cover of the disk with three punctures. The double cover branches at each of 
the punctures and at the disk's outer boundary. The interchange of rods given by cri and o<i 



2 The computer simulation were performed using Matthew D. Finn's code for solving Stokes' equation. This 
code uses complex variable methods to guarantee an accurate solution. 
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FIG. 3. The toral LTM consists of two linear shears. The first shear (left image) is in a vertical 
strip on the torus. The left edge of the strip remains fixed, while the right edge is translated by 
k units. The second shear (right image) is in a horizontal strip. The bottom edge remains fixed, 
while the top edge is translated by £ units. Here, the shears are shown with k = £ = 1. 



in the three-rod system correspond to the vertical and horizontal Dehn twists in the toral 
LTM. 

The domain of the toral LTM is an L-shaped subset of the flat torus formed by two 
overlapping orthogonal strips. One full application of the map consists of two consecutive 
linear shears. The first shear takes place in the vertical strip of width a, and the second 
shear is in the horizontal strip of height (3 (Figure [3]). The integer parameters k and £ are 
the strengths of each shear. In the limit a — (3 — 1, toral LTMs reduce to toral maps with 
uniform stretching. These are often called generalized Arnold cat maps: the cat map is 
the toral LTM with a = (5 = k = £ = l. Toral LTMs can also be defined more generally 
with non-linear shears [18J, but we do not deal with those here. If k£ > we say that the 
system is counter-rotating; if k£ < we say that it is co-rotating. This naming convention 
is the opposite of the one used by some other authors, for instance [IB]. In the present 
context, given the rod-stirring applications, it is more natural to define co-rotating and 
counter-rotating as we do here. 

Although we don't have a true fluid flow here, we can still compare the topological entropy 
measured from the stretching rate of a material line to a lower bound for the topological 
entropy of the map. This lower bound arises from a semi-conjugacy between the toral LTM 
and the generalized Arnold cat map with the same k and £ [27]. The generalized Arnold 
cat map acts on the full torus with the same matrix given in ([T]). For k and £ such that the 
matrix is hyperbolic, the topological entropy is again the log of the magnitude of the largest 
eigenvalue. 

As a parallel to our (J\(J^ and o\o\ examples above, we look at toral LTMs with k = 1, 
£ = 1 (counter-rotating) and k = 1, £ = —5 (co-rotating). Figure [I] shows, for both 
the counter-rotating and co-rotating maps, the image of an initial line segment after two 
iterations. As with the hydrodynamic flow examples, they both have a lower bound of 
^hom = 0.962. Numerically measuring the stretching of the line yields h^ ow = 0.962 for 
the counter-rotating LTM and ha ov = 1.91 for the co-rotating case. As for the three-rod 
mixers, the entropy agrees well with the lower bound for the counter-rotating example (< 1% 
difference), but there is a large gap for the co-rotating example (/ihom accounts for roughly 
49.6% of /iflow)- Looking again at Figure [ij we notice visible secondary folding in the co- 
rotating case (indicated by boxes), but none in the counter-rotating case. This holds in 
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(a) counter-rotating (b) co-rotating 



FIG. 4. Images of the same initial line segment after two iterations for (a) a counter-rotating toral 
LTM (k = £ = 1), and (b) a co-rotating toral LTM (k = 1, £ = 5). The small boxes highlight 
examples of secondary folding (acute angles) in the co-rotating case. 
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FIG. 5. The three regions given by Equation ([2]), shown here with k = 1. 
general for toral LTMs. To prove this, we take a closer look at the toral LTM map. 

A. Secondary Folding in Toral LTMs 

Recall that the toral LTM is the composition of a shear in a vertical strip of the torus 
and a shear in a horizontal strip on the torus. As a result, the domain is partitioned into 
three regions: points that undergo the vertical shear, but not the horizontal; points that 
undergo the horizontal shear, but not the vertical; and points that undergo both shears. We 
label these three regions Rv, Rh, and R^ as follows (see Figure O: 

Rv = {(x, y) mod 1 1 < x < a, (3 < y + xkoT x < 1}; (2a) 
R H = {( x , y) mod 1 1 a < x < 1, < y < /?}; (2b) 
R L = {( x , y) mod l|0<x<a, 0<y + xka' 1 < /?}, (2c) 

and also define R = Rv U Rh U -ftr • More specifically, for a point z = (x, y) T , the toral LTM 
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map is given by 



M(z) 



Vz mod 1, if z G Ry] 
Hz mod 1, if z G -R#; 
Lz mod 1, if z G Rl, 



(3) 



where V, H, and L are the matrices 



V 



1 

« 1 



H 



1 A 
1 



HV 



1 + kA A 
« 1 



(4) 



with k = fca -1 and A = The toral LTM is continuous and piecewise linear. 

Because of the piecewise nature of linear toral LTMs, the image of a line segment that 
crosses from one region to another will be a segmented line. This is the origin of all bends 
in iterates of line segments. Bends can be either obtuse or acute. We conjecture that obtuse 
bends do not contribute to the topological entropy, and so we ignore them; however, acute 
bends (also called kinks) always appear associated with an increased topological entropy. 
These kinks are the manifestation of secondary folds on the torus. To show the existence (or 
non-existence) of kinks, we examine the unstable manifold. This is sufficient since, under 
repeated iteration, a line segment will align with the unstable manifold and stretch in that 
direction. 

A subtle but important point here is that there is a set of singularities, S G R, where the 
unstable (and stable) manifolds don't exist [TBI EE] ■ This set consists of all points on the 
boundaries of the Rl, Ry, and Rh regions, as well as all (forward and backward) images 
of these points. These singularities exist because the toral LTM is continuous, but not 
differentiable, on the boundaries of the regions, so the tangent vector is not well-defined 
at those points. (These singular points eventually become the vertices of bends in line 
segments.) Fortunately, S is a set of measure zero [28], and tangent vectors are well defined 
in the neighborhood of the singularities. So, using continuity, we will speak of having an 
'unstable' manifold at each singular point, even though it may be the vertex of a bend. 
Then, since line segments tend to align with the unstable manifold, if a system has a kinked 
unstable manifold then any line segment will also develop kinks. 

The following lemma gives the slope of the unstable manifold at any (non-singular) point. 

Lemma 1. For a toral LTM, the slope of the unstable manifold at a point z G R\S is given 
by 

SJz) = — (5) 



Ani + 



1 



KlTLi + 



1 



\n 2 H ; 

Km 2 + ... 

where mi,rii > if z is in the horizontal strip, and mi,Tii > but rii — otherwise. The rrii 
and ni are unique integers that result from iterating z backwards along its orbit. 

(Note that this continued fraction was derived in a slightly different manner in [28J.) 

Proof. We use the fact that under repeated forward iterations, all line segments converge to 
the unstable manifold. To approximate the unstable manifold at point z, we first follow z 
backwards along its orbit to some point z'. Then, we can take (almost) any vector in the 
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tangent space, T z /R = M 2 , and map it forward to z. The image of the vector at z will be 
approximately tangent to the unstable manifold at z. (This fails only when the initial vector 
is exactly aligned with the stable manifold.) 

Each iteration of the inverse linked twist map, M~ 1 {z), leads multiplication by either V^ 1 , 
H" 1 , or L^ 1 = V^H" 1 , depending on the location of z. So repeated backward iterations 
can be written as 

z > = M~ N {z) = V~ m iH- n t . . . V- mi H- ni {z) (6) 

where N = Xli=i( m « + n « ~ !)• (The value of j depends on both N and z.) In general, 
rriijUi > 0. However, for some points, the first action of the inverse map is V -1 , not H~ x , 
so for these points, n\ = 0. These are exactly the points that are not in the horizontal strip. 
(Additionally, for a given N and z, we might have rrij = 0. In this case, ^ below ends 
with Xrij] this does not change the argument.) 

Let w' be a vector in T Z >R = IR 2 . Since DV = V and DH = H, w and z' evolve according 
to exactly the same matrix multiplication. That is, w G T Z R can be written as 

w = H ni V mi H n2 V m2 ... H n i V m >w' . (7) 

Now consider how multiplication by H n V m affects the slope of a vector. If the vector initially 
has slope s, then after multiplication by H n V m it has slope 

, s + Km 1 , . 

s = r = • 8) 

1 + Xn(s + Km) 1 v ' 

Xn H 

Km + s 

Repeated iteration to find the slope at w gives the finite continued fraction 



slope(w) = . (9) 

Ari! H 

ami H 



Xn 2 



nm 2 



1 

+ — 



nmj 



To find the exact slope, let iV — > oo. This means j — > oo and the finite continued 
fraction becomes the infinite continued fraction given in ([5]). So it only remains to show 
that the continued fraction converges. In the counter-rotating case, k,X > 0, and we have 
Km; > 1 and Xrii > 1. Hence ^](«;mj + Xrii) diverges, and the continued fraction converges 
to a finite value by the Seidel-Stern Theorem [29J. 

By definition, the co-rotating toral LTM must have k£ < 0, or equivalently kX < 0. 
However, we add another restriction and require kX < —4. This condition ensures that 
matrix L is hyperbolic and, consequently, that the toral LTM has an ergodic partition [28J . 
(With a slightly stronger condition, the toral LTM will also be Bernoulli [19J.) To show 
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convergence in this case, we transform the continued fraction into the form 

1 + - 



(KAm 2 n 2 ) 1 

1 + 



1 + 

Note that since k\ < —4, each of the numerators (except for the first) has magnitude less 
than 1/4. The continued fraction then converges to a finite value by Worpitzky's Theo- 
rem [29J. (The case where ni = follows similarly.) □ 

We intend to show that counter-rotating toral LTMs do not have kinks, while co-rotating 
toral LTMs do. Therefore, we proceed by investigating each type separately. For conve- 
nience, we assume in both cases that k > 0. 



1. Counter-rotating 

Since we are only considering maps with k > 0, we have £ > as well. Then every part 
of the continued fraction (jsl) is positive, and we conclude that the slope of the unstable 
manifold, S u , is positive at every point z G R\S. Furthermore, it is easy to show that the 
orientation of the unstable manifold is preserved under the action of the map. 

Define two cones in tangent space TR = M 2 by 

G\ = {(u,v) | u > 0,v > 0}, (11) 
C 3 = {(u,v) | u < 0,v < 0}. (12) 

(The subscripts refer to the standard quadrant numbering.) Note that the definitions are 
independent of z. 

Lemma 2. For k,£ > 0, the cones C\ and C% are invariant under both the vertical twist 
(V ) and the horizontal twist (H). 

Proof. Take a vector w = (u, v). Then Vw = (u, v + uk), and Hw = (u + vX, v). Clearly, if 
w G C\ then both Vw and Hw are also in C\ (and therefore hw G C\). Similarly, if w G C 3 
then both Vw and Hw are also in C3. □ 

Corollary 3. A line segment with positive slope will never kink under applications of M. 

Proof. Consider a line segment, 7, with positive slope passing through a point z. Let W\ 
and W3 be vectors in T Z R that are tangent to 7 and such that w\ G C\ and w% G C3. 
Then from Lemma |2j DM(wi) G C\ and DM(w%) G C3 regardless of which region z is in. 
Consequently, the angle between the vectors is obtuse, and no kink has formed. □ 

Theorem 4. The unstable manifold of a counter-rotating toral LTM has no kinks. 

Proof. Suppose that there is a kink in the unstable manifold. Then there must be a portion 
of the unstable manifold that is a straight line, crosses a region boundary, and maps to the 
kink. But S u > for counter-rotating toral LTMs, which contradicts Corollary |3j Hence, 
there are no kinks in the unstable manifold. □ 
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2. Co-rotating 

Next we consider co-rotating toral LTM, and again restrict attention to ones with k > 
and k\ < —4, so that the matrix L is hyperbolic. For such toral LTMs, we prove that the 
unstable manifold is kinked. 

We begin by using the eigenvectors of L to define two pairwise-invariant cones (as in [IB]). 
For each point z G R, we define two cones in the tangent space T Z R by 

C(z) = {(u, v)\m*< v/u < 0}, (13) 
C(z) = {(u,v) | < m* + k < v/u}, (14) 
where m* is the slope of the expanding eigenvector of L and for co-rotating maps: 

2k kX + ^kX(kX + 4) 



kX - a/kA(kA + 4) 



2A 



< 0. (15) 



Note that cone definitions are independent of z, so we can refer to them simply as C and 
C . For convenience, we let Iq = [to*, 0] and Ig = [m* + k, oo) denote the range of slopes in 
cones C and C respectively. The following lemma shows that tangent vectors to the unstable 
manifold lie precisely in these cones. 

Lemma 5. For kX < —4 and z G R\S, the slope of the unstable manifold, S u (z), lies in C 
if z is in the horizontal strip, and lies in C otherwise. 

Proof. This follows from a modification of a theorem by Hillam and Thron regarding conver- 
gence regions of continued fractions [30]. They look at a sequence of maps, Fi = /io/ 2 o- • -of i} 
where fi has the form fi(r) = ai/ibi + r), and the sequence {Fj(0)} gives successive conver- 
gents of a continued fraction. They prove that if there exists a disk D = {w G C | \w— c| < r}, 
with |c| < r and fi(D) C D, then the continued fraction converges to a value in D. 
We modify this proof for our purposes so that /,(r) has the form 

fi(r) = (16) 

Xni H 

nrrii + r 

Now the sequence {Fi(0)} gives the even numbered convergents of the continued fraction 
(j5|. Since we know that the continued fraction converges (Lemma [T]), these even numbered 
convergents will converge to the same value as the continued fraction. 

The disk we use is D = {w G C | \w — rn if /2\ < — m*/2}. This is the smallest disk 
possible, because when rij = rrii = 1, /i(m*) = m*. To see that fi(D) C D, re-write fi(r) as 

ArijT + nXniirii + 1 

and view it as a Mobius transformation. Then fi(dD) is a circle with center and radius 
given by: 

1 / KTUi + KfUi \ 

center = - 1 — (18) 

2 \1 + KXiriiUi 1 + Xni{m* + Krrii) J 

radius = — — — — -r. (19) 

2(1 + KXrriini) 1 + XnAm* + «mj 
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FIG. 6. A (possibly jointed) line segment passing through a point on the boundary of and Rh, 
with slope m 6 (— n — A -1 , 0], will form a kink after one iteration of M. 



From these, it is easy to check that fi(D) C D. The intersection of D with the real axis is 
the interval Iq = [m,*, 0], i.e. exactly the range of slopes in C. Thus, for rii > 0, S u (z) G Ic 
and the unstable manifold is in C. 

When ni = 0, equation ^ has the form 

S u (z) = Km, + — (20) 

Xn 2 H 

nm 2 + 



Xn 3 + ... 

and S u (z) G follows in the same manner. □ 

The following lemma shows how a kink may be formed after one application of the toral 
LTM. 

Lemma 6. Take a point z G S C R on the boundary between the and Rh regions (i.e. 
z G Rl H Rh)- Then a line segment passing through z (possibly with a bend at z) and with 
the segments on either side of z having slope s G (— K — A _1 ,0] will form an acute angle, 
with vertex at M(z), after one application of M. 

Proof. Suppose for simplicity that z is on the right-hand boundary of Rl, as in Figure [6} 
(The case where z is on the left-hand boundary is proved in an identical manner.) The 
portion of the line segment that lies to the left of point z is in the region Let u = (u±, u 2 ) 
be a vector that is parallel to the line segment and points left. That is, u\ < and, because 
u 2 /ui G (—ac — A -1 ,0], we have u 2 > 0. In this region the action of the toral LTM is 
multiplication by the matrix L = HV . So if the segment left of z is initially parallel 
to u, then it will be parallel to u' = Lu after M has been applied. We can compute 
u' = (ui + U\K,X + u 2 \,u 2 + Uik), and it is easy to see that u' points strictly into the 4th 
quadrant. 

Meanwhile, the portion of the line segment that lies to the right of point z is in the 
region R H . Let v = (vi, v 2 ) be a vector that is parallel to the line segment and points right. 
That is, V\ > and v 2 < 0. Since we are in region Rh, the action of the toral LTM is 
multiplication by the matrix H. So after M is applied, the line segment will be parallel to 
v' = Hv. Simple computation yields v' = (t>i + v 2 X, v 2 ), and we can see that v' points either 
horizontally or into the 4th quadrant. Thus, M(z) is the vertex of an acute angle and a kink 
has formed. □ 
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Theorem 7. The unstable manifold of a co-rotating toral LTM with k\ < —4 has kinks. 

Proof. Lemma [5] shows that for non-singular points in the horizontal strip, the unstable 
manifold is in cone C. Since Iq C (— k — A~ x , 0], pieces of the unstable manifold that cross 
the boundaries between R L and Rh satisfy the hypotheses of Lemma [6j Thus, these pieces 
kink after one application of M. Since pieces of unstable manifold map to other pieces of 
unstable manifold, this implies that there already existed kinks in the unstable manifold. □ 



IV. DISCUSSION 

We set out to explain the gap between the homological lower bound on the topological 
entropy and its numerically-observed value. We observed that, in both fluid-dynamical 
systems and toral LTMs, the presence of a gap appeared associated with 'secondary folding' 
- the presence of extra folds in material lines, not associated with topological obstacles such 
as rods. Though we have not been able to rigorously show that the gap is due to the folds, 
we were able to show that counter-rotating toral LTMs never have folds, while hyperbolic 
co-rotating toral LTMs always have folds. This correlates perfectly with the presence or 
absence of a gap in the topological entropy. Future work will aim to show that secondary 
folding is the direct cause of the extra entropy. 
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